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Abstract 

We investigate the zero-temperature quantum phase transition of the random 

bond Ising chain in a transverse magnetic field. Its critical properties are 

identical to those of the McCoy- Wu model, which is a classical Ising model 

in two dimensions with layered disorder. The latter is studied via Monte 

Carlo simulations and transfer matrix calculations and the critical exponents 

are determined with a finite-size scaling analysis. The magnetization and 

susceptibility obey conventional rather than activated scaling. We observe 

that the order parameter- and correlation function-probability distribution 

show a nontrivial scaling near the critical point which implies a hierarchy 

of critical exponents associated with the critical behavior of the generalized 

correlation lengths. 
75.10Hk„75.50Lk 
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Quite recently there has been a growing interest in the zero-temperature quantum critical 
behavior of disordered spin systems. Thermal fluctuations are absent here and the phase 
transition is driven by the interplay between randomness and quantum fluctuations. In 
order to tune the system to criticality one can either vary the strength of the disorder, 
as for instance in spin-| XXZ-chains JT], or one can control the strength of the quantum 
fluctuations directly by an external transverse magnetic field in spin models with a strong 
Ising anisotropy. The latter case is particularly interesting since experimentalists became 
able to investigated the effect of a transverse field on the glass transition in the Ising spin 
glass LiHo x Yi_ x F4 at low temperatures ||. On the theoretical side much progress has been 
made since then: the infinite range model has been solved analytically a Migdal- 

Kadanoff renormalization group calculation has been made 0, the critical exponents in 2 
and 3 dimensions have been determined via Monte Carlo simulations |7]|| and new results 
for the one- dimensional case have been derived via a renormalization-group analysis ||. 

In the latter papers focusing on the critical behavior in finite dimensions of these quantum 
models it has been pointed out that their universal properties are identical to those of clas- 
sical Ising models with layered disorder |TIJ. Especially the Ising spin chain in a transverse 
field can be mapped onto the McCoy- Wu model []TT]^T3[] , for which various exact results have 



been derived. The critical exponents of the order-parameter {3 and the correlation length v 
as well as the dynamical exponent z has been found recently via a renormalization-group 
(RNG) analysis by D. Fisher @. 

The aim of the present paper is to perform a numerical investigation of the finite-size 
scaling behavior of the Ising spin chain in a transverse magnetic field. Such an analysis has 
not been performed yet for this model and bears some new features concerning the finite- 



size scaling of anisotropic systems |T^]. It can also be seen as a test-ground for numerical 
methods applied to systems, for which - - in contrast to this model - - no quantitative 
theoretical prediction are at hand (e.g. for the cases considered in And finally it 

provides a check to what extend analytical predictions, like those made in |5| and which are 
valid asymptotically for an infinite system with rather unusual properties, can be detected 



in systems of finite size. 

The model under consideration is described by the quantum Hamiltonian 

^ = -E^m-rE<, (i) 

i i 

where a are spin— | Pauli matrices, T is the transverse field strength and the exchanges Jj are 
quenched random variables obeying a distribution P{J). At zero temperature the system 
(H) has a ferromagnetic phase transition to long-range magnetic order at a critical value 
r c , which depends on the bond distribution P{J)- We are interested in the critical proper- 
ties of this transition. However, as shown in P, |ll|,^3| , the magnetization behaves already 
non-analytically at higher values of T c giving rise to e.g. a divergence of the longitudinal 
susceptibility at higher values of T. 

The ground state energy of this one-dimensional quantum model ([!]) is equal to the free 



energy of the two-dimensional classical Ising model [15 



H — — ^ JiSijSi + \j — K ^ SijSij + i (2) 

at a certain finite temperature T. Here Sij = ±1 are classical Ising spins, the site index i 
runs along the x (space) direction and the index j along the r (imaginary time) direction of 
a two-dimensional square lattice. Following [|7],U we can rescale the bond strengths Jj and 
coupling constant K without changing the universal properties. For numerical convenience 
we set K = 1 and take a binary distribution 

P(J) = l6(J-j 1 ) + ±6(J-j 2 ), (3) 

where we put j\ = 1. The layered random bond Ising model (0) with the bond distribution 
(0) has a ferromagnetic phase-transition at a critical temperature T c defined by [[IT 



log coth (1/T C ) + log coth (j 2 /T c ) = 4/T c (4) 

and the universal properties (like exponents etc.) are identical to those of the quantum chain 
(H) at T c and zero temperature. Therefore we study model (0) at T c (j2) by Monte-Carlo 



simulations of rectangular lattices of size L x L T . The largest size in the r-direction was 
L T = 160 whereas that in the space direction x, which corresponds to the length of the 
quantum chain ([[]), was L = 16. Hence the disorder average over the distribution (||) could 
be done exactly by generating all, non-equivalent bond-configurations (whose number is 
approximately 2 L /L). To have more confidence on the data from the Monte-Carlo simulation 
we have compared our results to those obtained from transfer matrix calculations. The 
advantage of the latter method is that the results are exact, but the drawback is that one is 
limited to small system sizes L < 10. In all cases we found no significant deviations between 
the results of the two methods. 

The correlation length of the quantum chain (|l|) diverges as £ ~ (T — T c )~ u when ap- 
proaching the ferromagnetic transition. The characteristic relaxation time of the quan- 
tum dynamics is expected to diverge as r ~ £ z , with z being the dynamical exponent. 
These two diverging scales can naturally be found in the classical model (H): due to the ex- 
treme anisotropy one expects the correlation length in the space (or x) direction to diverge 
like £ ~ (T — T c )~ u and the correlation length in the imaginary time (or r) direction like 
£ T ~ (T— T c )~ zu . Following a nice argument by L. Mikheev |T^| one might imagine the system 



close to T c as being composed of roughly rectangular, ferromagnetic/paramagnetic domains, 
which are located (in the space direction) at segments of the chain where strong/weak bonds 
are dominating. The correlation length £ r in the time direction is then given by the average 
distance of domain walls in semi-infinite strips of width £ being ferromagnetically ordered, 
thus £ T oc exp(a£) (see e.g. |17]])- This supports very much an activated dynamics scenario 
with z = oo, as found in the RNG analysis 0. However, we begin the analysis by assuming 
a finite z and will see how far we get. 

At the critical point (i.e. T = T c ) various thermodynamic quantities are expected Jl4j to 
depend only on the scaling variable L/L T , the aspect ratio or the shape of the system. For 
instance we would have for the averaged spontaneous magnetization 

M = [<m)U « L-^m(L T /L z ) , (5) 
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where m = (L T L)~ l \Y2ij Sij\ is the magnetization per site, (• • •) means the thermal average 
and [■ ■ -] av means the disorder average. In figure 1 we show a scaling plot according to ([Sp 
obtained by Monte-Carlo simulations for different shapes and sizes (note that the disorder 
average is done exactly) at ji = 0.1 and T = T c = 1.32038. It yields (3/v — 0.17 ± 0.01 
and z = 1.65 ± 0.05. We also looked at j% = 0.05 for which T c = 1.14710, and obtain 
z = 1.70 ± 0.05 and (3/u = 0.18 ± 0.01. The RNG prediction is u = 2 (for the averaged 
correlation length) (3 = (3 — v^5) ~ 0.38 fU, which yields a ratio (3/v « 0.19 that agrees 
roughly with our estimate. 

The exponent z decreases as j% increases approaching z = 1 for j% = 1 and indicating 
the crossover to the pure case. This is the reason why we used small values for j 2 to ensure 
that what we see is the critical behavior of the disordered model. On the other hand we 
tried to avoid too small values of j2 in the Monte-Carlo simulations, since then the critical 
temperature decreases too much and equilibration becomes more difficult. 

The insert of figure 1 shows the scaling behavior of the susceptibility at j 2 = 0.1, 

X (L, L T ) = L T L[(m 2 )] av « L^'^x{L T /L z ) , (6) 

for which we find j'/u = 2.3 ± 0.1 (we used the value z = 1.65 that was obtained from the 
scaling of the magnetization). Note that 7' is not the critical exponent that describes the 
divergence of the susceptibility in an infinite system by approaching the temperature T c from 
above, since this quantity is expected to diverge already at a temperature higher than T c 



PJTT[| . The prediction for an anisotropic system obeying hyperscaling JTJj] is 7'/ v + 2(3 /v = 
d + z. Inserting d = 1 and the values for (3/u and z given above this relation is fulfilled very 
well. For the magnetization and the susceptibility conventional scaling seems to work well 
for these system sizes. 

The averaged cumulant g av = 0.5 • [3 — (m 4 )/ (m 2 ) 2 ] av is expected to scale like 

g av (L,L T )=g(L T /L z ) (7) 

and a scaling plot is shown in figure 2 with z = 1.55 ± 0.05 for 22 = 0.1, which is slightly 
lower than the estimate from the spontaneous magnetization. Furthermore the data collapse 



is not as good as in figure 1. Even worse is the scaling behavior of the cumulant g — 0.5 • 
(3 — [(^) 4 ]av/[(^) 2 ]av)- A systematic shift in the maximum to smaller values for increasing 
system sizes L indicates that this quantity is not dimensionless as expected naively. The 
natural scaling assumption [(m) fc ] av ~ L~ k ^l v m^[L T j L z ) does not seem to be correct here. 
As a consequence, the order parameter probability distribution P(M) = [S((m) — -M)] av 
does not scale in a trivial way like P(Ai) ~ L^^ U P{M.L^^^L T /L Z ) — as it does e.g. in 



conventional spin glasses [18j — which we checked explicitly by looking at the magnetization 
histograms for systems with constant aspect ratios L T m L z . 

To have an independent check of this scenario we have compared the above results with 
those obtained from the transfer matrix calculation. We have used the method recently 



introduced in |T^] for the exact calculation of the free energy derivatives. The method 
has been extended to a finite rectangular lattice with periodic boundary condition in both 
directions. By expressing the cumulants of the magnetization as derivatives of the free 
energy, we have computed g av and ~g for systems sizes up to 8 x 256. Since there are no 
numerical derivatives involved and the average over disorder is performed by summing over 
all possible configurations, the calculation yields the exact values. In all cases we did not 
find significant difference between these results and those of the Monte-Carlo simulations. 
We have also calculated the averaged spin-correlation function at T c , which is defined as 

C(r,t) = [(iSjjiSi +r j + £)] av . (8) 



For the averaged correlations in the time direction C(0,t) one expects |14]] for L T oc L z a 
behavior 

C(0, t) oc + (L T - t)r ni - , (9) 

where the second term on the r.h.s. takes into account the periodic boundary conditions. In 
the insert of figure 3 we have depicted C(0, t) for various system sizes L with L T chosen at 
the maximum of (7 av (L T ), so that L T oc L z . From the fit we conclude that 7]± = 0.23 ± 0.01. 
Concerning the spatial correlation function C(r, 0) Shankar and Murthy report a result 



[see equations 3.39 and 3.43 in |13|), which is 



C(r, 0) oc exp(— ay/r) + exp(— a\J L — r) , (10) 

where the second term on the r.h.s. takes again into account the periodic boundary condi- 
tions. This form yields a nice least square fit to our numerical data, as shown in figure 3. 
Note that comparing (||) and (|T0|) one observes that "space" and "time" seem to scale like 
r ~ (logt) 2 , as predicted in |J. 

D. Fisher argues that the result (|10D holds for the typical correlations (i.e. those 
corresponding to the maximum of their probability distribution), whereas the average should 
decay algebraically, since the latter should be dominated by rare, strongly correlated regions 
of the spin chain. We tried to fit C(r, 0) to an algebraic decay similar to (|^) with a different 
exponent m, which gave much worse results. However, stipulating 

C(r, 0) « r~^c(r/L , L T /L Z ) , (11) 

the relation C(L/2,0) ~ L~ v ^c(L T / L z ) should hold, which gives indeed an acceptable data- 
collapse for r/|| = 0.40±0.02 (using our estimate z = 1.70). The latter result agrees well with 
the result 7/|| « 0.38 obtained in 0. Furthermore it is consistent with the scaling relation 
V±/V\\ = z (within the errorbars) and 7711 = d + z — "f/u = 2(3 jv. We would like to stress 
that our data are compatible with both equations ( |IDD and (|TTj), which, however, are based 
on assumptions excluding each other. 

Next we study the probability distribution of correlation functions in the spatial direction 
through the analysis of the associated generalized correlation lengths [p0| , pl]1 . The analysis is 
carried out using the transfer matrix approach, for more details see [pljp2[ . We focus on the 
spatial correlation function for which we can use the general results for products of random 
matrices. In this case the latter consists of a succession of transfer-matrices from row % to 
row i + 1 in the spatial direction, each of which is made by L T spins. Since the generalized 
correlation lengths to be defined below are related to the first two Lyapunov exponents of an 
infinite product of transfer matrices we have to choose L very large (~ 10 6 ). An advantage 
of this approach is that one is left with only one scaling variable at the critical point since 
L T /L Z is zero. 
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In general each moment of the probability distribution of the correlation function defines 
a characteristic length scale which we denote by £ 9 where q is the order of the moment 

IB 



C 1 = - lim — lnfGifr) 9 ! (12) 

where Gi(r) is the connected correlation function between the row i and the row % + r. 
For example £1 is the characteristic length scale of the average correlation function, while 
£o is that of the typical correlation function. We shall then call £o the typical correlation 
length and £1 the average correlation length. It can be shown that if q > q' then £ g > £ g / 
|22|]. Taking into account this hierarchy, the usual finite-size scaling hypothesis jl4j for the 



spatial generalized correlation lengths would be 

C q l = L^^UiT -T c )L l J z ^) (13) 

where the z q can be called "generalized dynamical exponents". From the relation £ g > £ q i 
for q > q' it follows that z q > z q > and v q > v q i. In the scaling form ( |13"D it is assumed 
that the £ q diverges all at the same temperature T c , which in general need not to be the 
case, see for instance p0| . As consequence, especially for large q ( |13"D should be modified 
allowing for a g-dependent T c . However, we stick to the scaling form (pT3|) with T c given by 
the ferromagnetic phase transition temperature. 

We have calculated (for technical details see |21|J22| ) the z q by computing the exponents 



£g at the critical point T c for different ji for systems of sizes up to L T = 8. The length of 
the product was 10 6 . The exponent z q increases as j 2 decreases. We get (with at least 10 
percent accuracy) 
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In all cases we found a good scaling at the critical temperature T c given by (|). Our statistics 
is not accurate enough to investigate higher correlation lengths. The dynamical exponents 
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increase systematically with decreasing ji and it cannot be ruled out that z q — > oo for 
j 2 — > 0. Since the other exponents (as u q , see below, or (3 see above) are less susceptible to 
a variation in j'2 we have to leave it open here, whether this behavior indicates a cross-over 
or an actual non-universality of z q with respect to j%. 

The exponent v q are obtained from the data collapse using (|I3|). In this case we computed 
the generalized correlation lengths for temperatures T > T c and for system sizes up to 
L T = 7. Again the length of the product was 10 6 . With this statistics we were able to 
estimate only the first two exponents. The data for ji = 0.1, 0.05 and 0.01 leads to vq ~ 0.7 



and v\ ~ 1, in agreement with the analytical result of |T3[ but disagrees with the RNG result 
^=2|. In figure 4 we show a scaling plot for £ - 

We briefly discuss the time-correlation length for which we consider a semi-infinte strip 
in the time-direction [L T — > 00). While for the spatial correlations we have an infinite 
product of random matrices, in the case of time-correlations the transfer matrix is always 
the same since the randomness is only in the spatial direction. In this situation for each 
realization of disorder the correlation length is given by the inverse of difference between 
the first two eigenvalues of the transfer matrix and yields, if averaged, the inverse of the 
typical correlation length, ^q" 1 (note that this is equivalent to averaging the logarithm of the 
correlation function). For the typical time-correlation length we stipulate again the usual 



finite-size scaling form [14 



C^o = L~~ z ° ( T ,o((T — T c ) L 1 ^ ). (15) 

For j 2 = 0.1 we found, for example, zq ~ 1.3 and V m 1. Note that z 7^ z and v 7^ v. 

To summarize we performed a detailed finite-size scaling analysis of the zero-temperature 
phase transition occurring in a random bond Ising chain by tuning the transverse magnetic 
field to some critical value. For this model many analytical results are known and our 
analysis shows a good agreement with the results of Shankar and Murthy ||13| and concurs 
also with the RNG prediction of the existence of different length scales with different critical 
exponents v ||. However, we do not find the same values as those reported in the latter. 
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One possible explanation of this fact might be the following: D. Fisher || estimated the 
average correlations by only taking into account the very rare events, which, in our notation 
(12) for the generalized correlation length £ q corresponds to the limit q —>■ oo. Consequently 
his result v — 2 should an upper bound of our v q . What is denoted by typical correlations 
in |5| seems to us more related to our averaged correlation functions. 

Our finite-size analysis leads to finite values for the dynamical exponent z, which might 
be due to the small system sizes we were confined to. However, for j2 — > we find z — > oo in 
agreement with a RNG-picture ||. Furthermore we find that the order-parameter probabil- 
ity distribution scales non-trivially at the critical temperatures. We made this observation 
also in connection with the cumulants of the probability distribution of correlation functions 
and we found a hierarchy of critical exponents for the generalized correlation lengths. Despite 
these facts the numerical value (3/v for the finite-size scaling of the averaged spontaneous 
magnetization concurs with the prediction made in ||. 

All these phenomena merit further investigation (a more detailed discussion on these 
results will be given elsewhere [23|), especially with respect to the Griffiths singularities 
occurring already at temperatures above T c P,|T3|]. Finally we would like to mention that 
there is also a large overlap of the scenario we have encountered here with what might occur 
in two- and three-dimensional Ising spin glass in a transverse field . 



ACKNOWLEDGMENTS 

AC is grateful to G. Parisi, E. Marinari and A. Vulpiani for an interesting discussion and 
thanks the Institut fur Theoretische Physik of the University of Koln for hospitality where 
part of this work was done. HR would like to thank L. Mikheev and A. P. Young for many 
extremely valuable discussions. He acknowledges the allocation of about 10 5 processor hours 
on the transputer cluster Parsytec-GCell024 from the Center of Parallel Computing (ZPR) 
in Koln. This work was performed within the SFB 341 Koln-Aachen-Julich. 



10 



REFERENCES 



[1] C. A. Doty and D. S. Fisher, Phys. Rev. B 45, 2167 (1992); S. Haas, J. Riera and 
E. Dagotto, Phys. Rev. B 48, 13174 (1994); K. J. Runge and G. T. Zimanyi, |con^ 
mat/ 93 120021 

[2] W. Wu, B. Ellmann, T. F. Rosenbaum, G. Aeppli and D. H. Reich, Phys. Rev. Lett. 
67, 2076 (1991); W. Wu, D. Bitko, T. F. Rosenbaum and G. Aeppli; Phys. Rev. Lett. 
71, 1919 (1993). 

[3] J. Miller and D. Huse, Phys. Rev. Lett. 70, 3147 (1993); J. Ye, S. Sachdev and N. Read, 
Phys. Rev. Lett. 70, 4011 (1993). 

[4] For earlier work on the mean-field theory if quantum spin glasses see: Y. Y. Goldschmidt 
and P. Y. Lai, Phys. Rev. Lett. 64, 2567 (1990) and references therein. 

[5] For the mean-field theory of intinerant spin glass models see R. Oppermann, Nucl. Phys. 
B 401, 507 (1993). 

[6] B. Boechat, R. R. dos Santos and M. A. Continentino, Phys. Rev. B 49, 6404 (1994); 
M. A. Continentino, B. Boechat and R. R. dos Santos, preprint (1994). 

[7] H. Rieger and A. P. Young, submitted to Phys. Rev. Lett. (1994). 

[8] M. Guo, R. N. Bhatt and D. A. Huse, submitted to Phys. Rev. Lett. (1994). 

[9] D. S. Fisher, Phys. Rev. Lett. 69, 534 (1992). 

[10] For recent work on layered Ising models see L. V. Mikheev and M. E. Fisher, Phys. Rev. 
B 49, 378 (1994) and references therein. 

[11] B. M. McCoy and T. T. Wu, Phys. Rev. 176, 631 (1968); 188, 982 (1969); 188, 1014 
(1969) 

[12] P. Hoever, W. F. Wolff and J. Zittartz, Z. Phys. B 41, 43 (1981). 

11 



[13] R. Shankar and G. Murthy, Phys. Rev. B 36, 536 (1987). 
[14] K. Binder and J. S. Wang, J. Stat. Phys. 55, 87 (1989). 
[15] M. Suzuki, Progr. Theor. Phys. 56, 1454 (1976). 

[16] We would like to thank Lev Mikheev for a very enlightening discussion on this point. 
See also L. Mikheev, in preparation (1994). 

[17] P. Kleban and G. Akinci, Phys. Rev. B 28, 1466 (1983). 

[18] R. N. Bhatt and A. P. Young, Phys. Rev. B 37, 5606 (1988). 

[19] A. Crisanti, G. Paladin, M. Serva and A. Vulpiani, Phys. Rev. E, in press (1994). 

[20] B. Derrida and H. J. Hilhorst J. Phys. C 14, L544 (1981). 

[21] A. Crisanti, S. Nicolis, G. Paladin and A. Vulpiani, J. Phys. A 23, 3083 (1990). 

[22] A. Crisanti, G. Paladin and A. Vulpiani, Products of random matrices, Springer series 
in solid state sciences 104; Springer, Berlin-Heidelberg-New York (1993). 

[23] A. Crisanti and H. Rieger, in preparation 



12 



FIGURES 

FIG. 1. Scaling plot of the magnetization M(L,L T ) for j 2 = 0.1 at T c . It yields z = 1.65 ±0.05 
and j3/v = 0.17 ± 0.01 The insert shows the scaling plot of the susceptibility %(L, L r ) with 
ijv = 2.3 ± 0.1 and z as for M. 

FIG. 2. Scaling plot of the averaged cumulant g av (L,L T ) for 22 = 0.1 at T = T c with 
z = 1.55 ±0.05. 

FIG. 3. Insert: Correlation function in the (imaginary) time direction C(0,t) for 22 = 0.05 at 
T = T c . The system sizes are 4 x 16, 6 x 32, 8 x 64, 12 x 100 and 16 x 160, i.e. their aspect ratio is 
close to the maximum of the cumulant <7 av and therefore roughly constant. The full curve is a least 
square fit to C(0,t) oc + (160 - t)~^} and yields rf ± = 0.23 ± 0.01. Left part: Correlation 
function in the space direction C(r, 0) for 22 = 0.05 at T = T c . The system size is 16 x 160, close to 
the maximum of g av . The full line is a least square fit to C(r, 0) oc exp(— ar 1//2 )±exp(— a (L— r) 1 / 2 ) 
with a = 0.62. 

FIG. 4. Scaling plot for the typical correlation length in the spatial direction ^L^* versus 
(T — T c )L l J ZoUo . The upper data-set is for 22 = 0.05, where zq = 2.12 and vq = 0.7 has been used, 
the lower data-set is for 22 = 0.1, where zq = 1.74 and uq = 0.7. 
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